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{iLiI We investigate effects of threshold resummation of logarithmic correc- 



tions InA^ in Mellin space quantitatively. Threshold resummation leads to 
enhancement of next-to-leading-order QCD predictions for jet production at 
large jet transverse energy, which is in the trend indicated by experimental 
data. We show that this enhancement is completely determined by the behav- 
ed I ior of threshold resummation at small N, the region where hierachy among 
different powers of In A^ is lost and current next-to-leading-logarithm resum- 
mation is not reliable. Our analysis indicates that more accurate threshold 
resummation formalism should be developed in order to obtain convincing 
predictions. 
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1 INTRODUCTION 

The formalism of threshold resummation of double logarithmic corrections 
to QCD processes, which occur in extreme kinematic conditions, has been 
developed for some time |l|, 0, ||. At the kinematic end points, a special type 
of corrections ln(l— a;)/(l— x)+ is produced with x being a parton momentum 
fraction, which appears as ln^(l/iV) under the Mellin transformation. There 
have been abundant formal derivations of threshold resummation for various 
processes, such as deep inelastic scattering, Drell-Yan, direct photon, and 
heavy-quark productions [Q. Quantitative studies of threshold resummation 
effects in heavy quark production |^ and in direct photon production [Q 
have been performed recently. Such numerical studies are essential in order 
to justify that threshold resummation indeed collects important dynamics of 
processes in extreme kinematic conditions. 

In this letter we shall analyze effects of threshold resummation from an- 
other point of view. For our purpose, it suffices to consider resummation 
of the logarithmic corrections that can be factorized into parton distribu- 
tion functions (PDFs). Taking jet production at Tevatron as an example 
[0], we observe that threshold resummation enhances next-to- leading-order 
(NLO) QCD predictions for jet production at large jet transverse energy Ex. 
This tendency is qualitatively consistent with experimental data and with 
the conclusion in [^ |^. However, the enhanced predictions overshoot data 
by a factor of 2. A simple investigation reveals that the resummation associ- 
ated with the quark distribution functions lead to a negligible effect and the 
overestimation is attributed to the resummation associated with the gluon 
distribution function. The reason is that the color factor A''^ = 3 in the latter 
case is larger than Cp = 4/3 in the former case. 

We further find that the enhancement is completely determined by the 
behavior of the resummation associated with the gluon distribution function 
at low A^. Unfortunately, this is a region where current next-to-leading- 
logarithm (NLL) resummation is not reliable, since hierachy among different 
powers of InA^ is lost and all nonleading logarithms need to be summed. 
In other words, the large- A^ behavior of threshold resummation is reliable, 
but almost irrelevant to the end-point enhancement, whereas the small- A^ 
behavior accounts for the end-point enhancement, but is not reliable in NLL 
threshold resummation. Moreover, the importance of low-A^ contributions 
can not be diminished no matter how extreme kinematic conditions are. Our 



analysis indicates that more accurate threshold resummation formalism need 
to be developed in order to obtain convincing predictions. 

2 THRESHOLD RESUMMATION 



Threshold resummation of the logarithmic corrections that can be factorized 
into PDFs is written as [Hi 



/(iV) = exp 



1 , l-Z^-^ /•! d\ 

dz / -— 

1 — z J(i-z)2 A 



7x(as(v^p+)) 



(1) 



where p"*" is the longitudinal component of hadron momentum p, and the 
anomalous dimension '^k is given, up to two loops, by [^ 



7x = -C^+(-) Cp 

IT \ TV 



1 36 12/ 18 ^ 



(2) 



with Ca = 3 being a color factor and Uf = 4 the number of quark flavors. In 
this work we shall adopt its modifled version. 



f{N) = exp 



1-1/^ dz /-i 



1 — z 7(1-^)2 A 



--7i<:(a^(yAp+)) 



(3) 



which is equivalent to Eq. (|T]) up to 0{1/N) corrections. We refer readers 
to for its detailed derivation. Equation (H) is simpler in analytical ma- 
nipulation, since the integral in the exponent can be worked out explicitly. 
Threshold resummation associated with the gluon distribution function is 
obtained by substituting Nc = 3 for Cp in Eq. (|[). 

The first term of jk leads to the leading (double)-logarithm summation, 
and the second term leads to the NLL summation. Note that Eq. (^ is not 
complete at the NLL level. To obtain a full NLL summation, contributions 
from soft gluon exchanges among initial- and final-state partons should be 
taken into account, which are pro cess- dependent. While Eq. (|^) is process- 
independent, since it sums factorizable corrections. In this work we shall 
comment on applications of threshold resummation to various QCD pro- 
cesses, such as deep inelastic scattering, direct photon production, and jet 
production. Furthermore, our goal is to demonstrate that the end-point en- 
hancement is determined by the small- A^ behavior of threshold resummation. 
Hence, Eq. (H) serves the purpose. 



We propose an expansion of the threshold resummation f{N) in Eq. (|^) 
in terms of polynomials in N: 

fiN) = J2a.CiN-l,z), CiN,t)^--—-. (4) 

The above series corresponds to a simple form in momentum fraction space, 

f{x) = Y.^6^Hl-^), (5) 

1=0 

which can be easily verified by performing the Mellin transformation 

/(iV)= f'dxx^-'fix). (6) 







The first coefficient oq = /(I) = 1 gives the initial condition of threshold 
resummation. The other coefficients Oj are determined by best fit to Eq. (j^). 
On the other hand, Eq. (^ is not appropriate for large N, the region in which 
the integration variable A may be as small as l/N"^, the running coupling 
constant asip'^ /N) diverges, and perturbation theory is not applicable. In 
this region Eq. (^) should be replaced by a nonperturbative function, which 
is of course model-dependent. For example, a minimal prescription which 
takes into account only the Landau singularity was introduced in |T0|. It 



is straightforward to extrapolate Eq. (^) to the N -^ oo limit, and this 
extrapolation can be regarded as a nonperturbative model. 

It turns out that an expansion up to n = 3 in Eq. (P describes the growth 
of Eq. (^ with A^ very precisely. The parameters a, for the quark and gluon 
distribution functions from best fit to Eq. (|) for Aqcd = 0.2 GeV in a^ and 
for center-of-mass energy y^ = \/2p^ = 1800 GeV at Tevatron are listed 
below: 

parton ai 02 0.3 

quark 1.6198 x lO^^ -8.5872 x 10"^ 2.2515 x 10"^ 
gluon 1.1248 X lO^^ 7.2253 x 10"^ 2.5192 x 10"^ 

For smaller -^i, such as ^/s = 38.7 GeV for direct photon production in E706 
1^, we obtain the parameters 



parton ai 02 03 

quark 5.6075 x 10"^ 6.6109 x 10"^ -2.1474 x lO"*^ 
gluon 1.3616 X 10'^ 9.3512 x 10'^ 5.7660 x 10'^ 

The parameters Oj increase as ^/s decreases, implying stronger resummation 
effects, because the running a^ is larger at lower energies. For even lower 
^/s, such as those for deep inelastic scattering, Oj are even larger. 

The modified PDF (p is written as the convolution of threshold resumma- 
tion with the original distribution function (p |^: 

m = rf/(O0(^/o, 

Jx I; 

= (1 - ai + 02 - a3)0(x) - (ai - 2a2 + 3a3)x-—(f){x) 

ax 

1 d^ 1 d^ 

+-(a2 - 3as)x^—^ix) - g«3X^^0(x) . (7) 

Briefly speaking, threshold resummation effectively modifles a PDF, and the 
modiflcation is energy- and pro cess- dependent. Hence, before performing 
a global determination of PDFs, one should clarify threshold resummation 
effects. 

The motivation to expand the threshold resummation into a series of 
C{N-1, i)uptoi = 3 is as follows. The C{N-1, 1) term and the C{N-1, 3) 
term determine the behavior of the resummation at small N and at large N, 
respectively. If the series terminates at i < 3, the role of each C{N — l,i) 
in determining the behavior of the resummation in different regions of N is 
not signiflcant. If the series contains terms with i > 3, numerical handling of 
higher derivatives of modifled PDFs will be difficult. Since the resummation 
associated with the gluon distribution function dominates, we take it as an 
example to demonstrate the above idea. For A^ ~ 10, the i = 1 term aiC{N — 
1, 1) in the case with y/s = 1800 GeV is of order unity, while the i = 3 term 
a^C{N — 1, 3) is only of order 10~^. As A^ increases up to 10^, the edge 
for perturbation theory to be applicable, the i = 1 term, being of order 10^, 
becomes smaller than the i = 3 term, which is of order 10^. Hence, a variation 
of ai implies a variation of the small- A^ behavior of threshold resummation, 
and a variation of 03 implies a variation of the large- A^ behavior. 

For A^ ~ 10^, a double logarithm In A^ is larger than a single logarithm 
In A^ by a factor of 7, indicating that hierachy among different powers of In A^ 



exist and current NLL resummation is reliable. For N as small as 10, In^ N 
and In N are in fact of the same order, and NLL resummation is not reliable. 
Therefore, by varying the coefficients ai and 03, we can investigate how 
sensitive the end-point enhancement of NLO predictions is to the small-iV 
portion of the NLL resummation, which is not reliable, and to the large- 
N portion, which is reliable. It will be satisfactory, if the enhancement is 
insensitive to the small- A^ behavior of the NLL resummation. However, we 
shall demonstrate that this is not the case. The reason is obvious from 
Eq. (^: the coefficient of each term on the right-hand side of Eq. (0) is 
dominated by a, with smaller i. The above reasoning applies to cases like 
E706 with ^/s = 38.7 GeV or lower, for which the corresponding parameters 
Qi have the same relation ai ^ 02 ^ 03. In these cases the end-point 
enhancement is determined by the behavior of threshold resummation at 
smaller A^, and the controversy stated above is more serious. 

3 NUMERICAL ANALYSIS 

Consider jet production with transverse energy Et in pp collision, 

p{Pi)+P{P2)^AEt)+X. (8) 

The hadron momenta are assigned as pi = {pi,0,Ot) and p2 = (0,^2^, Oy) 
with Pi = P2 = V"^/^- Partons carry the momenta $,iPi with ^j, i = 1, 
2, being the momentum fractions. In fact, the transverse momenta kj^ of 
partons should be taken into account, when transverse degrees of freedom 
of final states are measured [0. The corresponding kx resummation, if 
included, is expected to further enhance the cross section at high Et- In the 
present work kr resummation will not be considered, since we concentrate 
on effects of threshold resummation. 

The factorization of jet production is basically similar to that of direct 
photon production in |jl2[. The self-energy correction to a parton and the 
loop correction with a real gluon connecting the two partons from the same 
hadron, contain both collinear divergences from the loop momentum / par- 
allel to Pi and soft divergences from small /. Since soft divergences cancel 



between the above corrections |T3[, the remaining collinear divergences are 
absorbed into a PDF associated with the hadron i. They are the corrections 
which produce the logarithms ln(l/iV) that have been summed into Eq. (^ in 
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axial gauge. In the considered kinematic region for jet production, the other 
radiative corrections, because of the soft cancellation, are absorbed into a 
hard scattering amplitude H, which corresponds to a parton-level differen- 
tial cross section. 

The factorization formula for jet production with the threshold resum- 
mation for PDFs included are written as 



daiEj 
dEr 



I d^id^2^{^i, Et/2)^{^2, Et/2)H{^,, 6, s, Et/2) . (9) 



We have set the renormalization (factorization) scale /x of and H to the 
characteristic scale Et/2. The original PDFs evolve to £^7-/2 according 
to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equation WM, which sums 



another type of single logarithms InEx- In the following numerical analysis 
we employ the NLO QCD calculations for jet production derived in []T3| 
and the CTEQ4M set [|T6| for the orginal PDFs. We have also checked the 
CTEQ3M set ^7\ and found that our conclusion does not depend on the 
choice of PDF sets. 

In Fig. 1 we show the modification from the threshold resummation on 
NLO QCD predictions for jet production at Tevatron. It has been observed 
that there seems to be an excess of the CDF data at high Et compared to 
the NLO predictions with usual PDFs [|1^, 0, |1^, 0, whereas the DO data 



are in good agreement with the predictions. Note that the CDF and DO data 
do not conflict each other if considering the large systematic uncertainties 
(8 ~ 30%) in addition to the statistical ones. Hence, it is expected that the 
threshold resummation causes a small amount of enhancement of the NLO 
predictions at high Et- However, it is found that the enhancement is a factor 
of 2. We explore the source that is responsible for this huge overestimation. 
If turning off the threshold resummation associated with the gluon distribu- 
tion function, the end-point enhancement falls dramatically. If turning off 
the threshold resummation associated with the quark distribution function, 
the enhancement almost remains invariant. That is, the resummation associ- 
ated with the quark distribution function contributes only few percent of the 
full enhancement, and the overestimation is attributed to the resummation 
associated with the gluon distribution function as shown in Fig. 1. 

We then investigate which portion in N of the resummation associated 
with the gluon distribution function accounts for the end-point enhancement. 
As stated before, the parameters ai and 03 control the small- A^ and large- A^ 
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Figure 1: Enhancements from the full threshold resummation, from the re- 
summation for quark only, and from the resummation for gluon only. The 
NLO predictions and experimental data (with statistical errors only) are also 
shown. 



behaviors of threshold resummation, respectively. If setting ai to zero, which 
changes the small- A^ behavior of the resummation but leaves almost invariant 
the large- iV behavior, the resultant predictions coincide with those from the 
NLO calculations without including resummation effects as shown in Fig. 2. 
If setting as to zero, which changes the large- A^ behavior but leaves invariant 
the small-iV behavior, the enhancement remains the same. This investiga- 
tion indicates that the behavior of the resummation at small- iV determines 
the end-point enhancement. As argued before, this is the region where hi- 
erachy among different powers of InA^ is lost and more accurate formulas 
including summation of all nonleading logarithms are required. Before ex- 



o 



o 
Z 
o 






2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 
0.8 
0.6 
0.4 
0.2 





CDF jets I — — 

DO jets ^-^- 
No resummation 
Resummation-gluon 

ai=0 

a3=0 - - - - 



%'»Tn9T^'^%^^^~'''T 



— 1 — ^. 



_l I I L. 



_l I I L. 



50 100 150 200 250 300 350 400 450 500 

E^. (GeV) 



Figure 2: Enhancements from the threshold resummation for gluon, from 
the resummation with Oi set to zero, and from the resummation with 03 
set to zero (which appears to coincide with dotted hne for gluon). The 
NLO predictions and experimental data (with statistical errors only) are 
also shown. 



tracting reliable predictions from threshold resummation, this point must be 
taken into account. At last, we observe that the end-point enhancement can 
be adjusted by tuning ai, i.e., by varying the small- iV behavior of thresh- 
old resummation. Choosing ai = 0.02 in the resummation associated with 
the gluon distribution function, the resultant predictions are similar to those 
from applying the CTEQ4HJ PDFs 0, |2^ , and well describe the CDF and 
DO data simultaneously. 

We emphasize that the above controversy always exists no matter how 
large Et is reached. With higher Et, behaviors of PDFs at larger momentum 



fraction are probed. In this region we have x — *> 1 in the inverse MeUin 
transformation, 

fix) ^ — / dNx-^'fiN) , (10) 

with c being an arbitrary constant, and contributions from f{N) in the whole 
range of A^ are equally important. In cases with intermediate Et, for which 
X may not approach unity, contributions from the small- A^ region dominate. 
That is, no matter how large Et is reached, the small- A^ region always con- 
tributes. We have also analyzed the data of direct photon production in 
E706 1TT| and of deep inelastic scattering in BCDMS ^ and NMC ^ with 



lower ^/s using the above method, and arrived at the same conclusion: in the 
region where perturbation theory is applicable, small- A^ contributions always 
determine the end-point enhancement. 

4 CONCLUSION 

In this letter we have shown that the behavior of the threshold resummation 
associated with the gluon distribution function at small A^ determines the 
end-point enhancement of NLO predictions for jet production completely. A 
variation of the large- A^ behavior of the NLL resummation, which is reliable, 
does not affect the enhancement. However, in the small- A^ region hierachy 
among different powers of InA^ disappears and current NLL resummation 
is not reliable. This controversy also exists in processes with higher E^ or 
lower center-of-mass energies. Hence, to obtain convincing predictions, more 
accurate formalism for threshold resummation which sums all nonleading 
logarithms need to be developed. We emphasize that it is not the goal of this 
work to explain experimental data. Even if data can be explained by complete 
NLL threshold resummation, the controversy we have found remains. Before 
attempting to understand data using threshold resummation, our conclusion 
should be taken into account. 



This work was supported by the National Science Council of R.O.C. under 
Grant Nos. NSC88-2112-M-006-013 and NSC88-2811-M007-026. We also 
thank National Center for Theoretical Sciences in Taiwan for partial support. 



10 



References 

[1] G. Sterman, Phys. Lett. B179, 281 (1986); Nucl. Phys. B281, 310 

(1987). 

[2] S. Catani and L. Trentadue, Nucl. Phys. B327, 323 (1989); Nucl. Phys. 
B353, 183 (1991). 

[3] G.P. Korchemsky and G. Marchesini, Nucl. Phys. B406, 225 (1993); 
Phys. Lett. B313, 433 (1993). 

[4] N. Kidonakis, G. Oderda, and G. Sterman, Nucl. Phys. B525, 299 
(1998); E. Laenen, G. Oderda, and G. Sterman, Phys. Lett. B438, 173 
(1998). 

[5] N. Kidonakis and R. Vogt, Phys. Rev. D59, 074014 (1999). 

[6] S. Catani et ai, |hep-ph/9903436i 

[7] CDF collaboration, F. Bedeschi, talk at 1999 Hadron Collider Physics 
Conference, Bombay, January, 1999; DO collaboration, B. Abbott et ai, 
Phys. Rev. Lett. 82, 2451 (1999). 

J.C. Collins and D.E. Soper, Nucl. Phys. B193 (1981) 381. 



[9] H-n. Li, |hep-ph/981 13401 , to appear in Chin. J. Phys; |hep-ph/9812363 



to appear in Phys. Lett. B. 

[10] S. Catani et ai, Phys. Lett. B351, 555 (1996); Nucl. Phys. B478, 273 
(1996). 

[11] E706 collaboration, L. Apanasevich et ai, Phys. Rev. Lett. 81, 2642 
(1998). 

[12] H.L. Lai and H-n. Li, Phys. Rev. D58, 114020 (1998). 

[13] H-n. Li, Phys. Lett. B369 (1996) 137; Phys. Rev. D55 (1997) 105. 

[14] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); G. 
Aharelli and G. Parisi, Nucl. Phys. B126, 298 (1977); Yu.L. Dokshitzer, 
Sov. Phys. JETP 46, 641 (1977). 



11 



[15] S.D. Ellis, Z. Kunszt, D. Soper, Pliys. Rev. Lett. 64, 2121 (1990). 

[16] H.L. Lai et al, Pliys. Rev. D55, 1280 (1997). 

[17] H.L. Lai et al, Pliys. Rev. D51, 4763 (1995). 

[18] A.D. Martin, R.G. Roberts, W.J. Stirling, Phys. Lett. B354, 154 (1995); 
A.D. Martin et al, Eur. Phys. J. C4, 463 (1998). 

[19] H.L. Lai et al, |hep-pli/990328^ . 

[20] J. Huston et al, Phys. Rev. Lett. 77, 444 (1996). 

[21] BCDMS collaboration, A.C. Benvenuti et al, Phys. Lett. B223, 485 
(1989). 

[22] NMC collaboration, M. Arneodo et a/., Nucl. Phys. B483, 3 (1997). 



12 



